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Abstract 

A new anisotropic t-J model in one dimension is proposed which has long- 
range hopping and exchange. This t-J model is only partially solvable in 
contrast to known integrable models with long-range interaction. In the high- 
density limit the model reduces to the XXZ chain with the long-range ex- 
change. Some exact eigenfunctions are shown to be of Jastrow-type if certain 
conditions for an anisotropy parameter are satisfied. The ground state as well 
as the excitation spectrum for various cases of the anisotropy parameter and 
filling are derived numerically. It is found that the Jastrow-type wave function 
is an excellent trial function for any value of the anisotropy parameter. 
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There have been remarkable developments in understanding the family of one- 
dimensional models with interaction proportional to the inverse square of the distance. [l|-[/J 
The simple structure of the spectrum places the family as canonical models for realizing 
the Tomonaga-Luttinger liquid in one dimension. Much less is known for the case where an 
energy gap is present for the excitation spectrum. The simplest example is the XXZ chain 
with the Ising anisotropy. It has been shown by Haldane |§ that for particular values of 
the anisotropy parameter, the long-range XXZ chain has the Jastrow-type wave function as 
an eigenfunction with non-zero magnetization. However, the wave function of the ground 
state is not of Jastrow type. The presence of the gap makes the system different from the 
Tomonaga-Luttinger liquid. In view of importance of the anisotropy in real physical systems, 
the study on the anisotropic models is urgently required. 

The purpose of this Letter is to present and solve a new type of the t-J model which 
reduces to the long-range XXZ chain in the high-density limit. We can derive a part of 
the eigenstates by mapping the model to the multicomponent Sutherland model. |7j The 
ground-state eigenfunction within a certain range of the magnetization is shown to be of 
Jastrow type. We clarify the condition under which the Jastrow-type function constitutes an 
eigenstate. For cases where analytic solution is not available, we diagonalize the Hamiltonian 
numerically. The ground-state wave function so obtained turns out to be extremely close to 
but not identical with the Jastrow-type one. As a preliminary to understanding of the full 
spectrum, we examine the excitation spectrum of the long-range XXZ model numerically. 

Let us first review the anisotropic spin model (XXZ model) with long-range exchange. |3j 
The Hamiltonian is given by 

N 

■H X Jy(S?S? + S?S ! j + A.V/.S:). (1) 

with Jij = JD(xi — Xj)~ 2 , where D(xi — Xj) = (N/tt) |sin[7r(xj — Xj)/N] \ is the chord distance 
consistent with the periodic boundary condition on N lattice sites. The parameter A(> 1) 
represents the Ising anisotropy. Haldane || has shown that the model is partially soluble 
under some conditions on anisotropy and magnetization, and has obtained a Jastrow-type 
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eigenfunction and its eigenenergy. The model we present below is an extension of the XXZ 
model with inclusion of holes. 

We introduce an anisotropic t-J model as follows: 



qx qx | QV CV i A Q z Q z 
Oj Oj T" *~>j *Jj I I— *2^i 



(A 3 - 1) (fi'fnj- + TiiS?) - -A 4 n i n i 



P, (2) 



where "P is the projection operator to exclude the double occupation at each site. We take 
Uj — Jij — tD(xi — Xj)~ 2 , where D(xi — Xj) is defined above. The parameters Ai and A2 
obviously break the SU(2) invariance. The other parameters A3 and A4 are necessary to 
make the model analytically solvable. It may seem highly unrealistic that the model has 
different magnitudes for the hopping. This difference originates simply from the mathe- 
matical requirement. If we translate the spin degrees of freedom into the chain index in a 
double-chain model, however, the difference is physically acceptable. When Aj = 1 for all i, 
this model reduces to the long-range supersymmetric t-J model. |J 

Any state in the Hilbert space can be represented by the positions of holes and down 
spins. If we let M denote the number of down spins and Q denote the number of holes, then 

iSf ot is given by S^ ot = (N — Q)/2 — M. The wave function is represented by 

M Q 

\1>)= E ^({ x h{v})I[s~l[hl.\F), (3) 

{x},{y} i=l j=l 

where \F) denotes the fully polarized up-spin state, S~. is the spin- lowering operator at site 
Xi, and hljj creates a hole at site yj. The amplitude ip({x}; {y}) is symmetric in the positions 
{x} = {x\, • ■ ■ , xm} of down spins, and antisymmetric in the positions {y} = {y%, • • • , uq) 
of the holes. When {x} and {y} are specified, the positions {u} = {ui, • • • , un-m-q] of the 
up spins are determined automatically. 

To solve the eigenvalue problem TC\ip) = E\ip), we divide @ as follows: Ti. — 7i' + TC rest , 
where 

H rest =?^[U(A 2 + 2A 3 - 2) (Sf + Sf) - \mA 2 + 2A 3 - 2) 



~ J ii (A 4 - 2A 3 + 2)(m + rij) + ^(A 4 - 2A 3 + 2 



This part 7i res t is diagonal i.e., TC res t\ 
eigenvalue. Namely we obtain 



V. (4) 
-Srestl'0), and gives a constant term E rest for the 



rest 
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(N 2 - 1) [( A 2 + 2A 3 - 2) (N - 4M - 2Q) - ( A 4 - 2A 3 + 2) (N - 2Q)} . (5) 



Next, we consider 7i' in the Hilbert space with the reference state \F). It is represented by 

w = - E E <*Mi - A i E E UjSrhthls} 

+ E E -V?^ + (A a + 2A 3 - 2) E E ^ - - \ 

ie{x}je{u] i£{u}mu},&i v z/ v 

-2(A 3 - 1) E E ^-n*nj fe - \) - ~ 



■~(A 4 - 2A 3 + 2) E E J ^ - - !)■ 



(6) 



In order to make the model soluble we parametrize Aj's as Ai = A, A 2 = A(A + l)/2, 
A 3 — 1 = A(A — l)/4, and A 4 = A(3 — A)/2 in terms of a single parameter A. Then 
TL\ip) = E\ip) turns into the following equation: 

M+QN-l 

E E J Mvi, ••-,% + r, ■■■ , v M+Q ) + E A(A + P l3 )J^({v}) — (E — E rcst )^({v}). (7) 

i=l r=l i<j 



where we use the notation {v} = {x} © {y}, and 



rest 
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(N - 1)[A(A - 1)(JV - 2Q) - 2A 2 M]. 



(8) 



With the plain wave basis, a general wave function can be represented as 



M+Q 

v>(M)= n ^ /2 e^ 

i=l {fc} 



M+Q 



(9) 



i=l 



where F of means the Fourier expansion and Z{ = exp(i27Tf i/iV). We shall specify the 
range of the momenta ki{i = 1, ■ ■ ■ , M + Q) of two kinds of particles (down spins and holes) 
shortly. Substituting (§) into (|7|), we obtain 



A^ 2 N 2 + 2 

7l z t O 



f (10) 
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Here we have used the relation: 



N-l A2irkr/N 



^ sm\irr/N) 



(ID 



for < k < N. In order to derive (|T0|) , it is thus required that —N/2 < ki < N/2. 
Since the effect of umklapp scattering is not considered in ip, we are not dealing with the 
complete set of wave functions. The left-hand side of (|T0|) is nothing but the Hamiltonian 
of the multicomponent Sutherland model with both fermions and bosons present. Thus 
the results of the multicomponent Sutherland model can be used to find the solution of the 
anisotropic t-J model. 

A Jastrow-type wave function 



M 



M 



M Q 



^(M; {y}) = n ^ JslN n ^ Jh/N n ^ - ^) A+l n ^ - %y n n ^ - %o 

j=l i<j i<j %=\j=\ 



i=l 



(12) 



is an eigenfunction of (0) where Aj's are parametrized with A, under the following conditions: 
A=odd integer, | J s — J h \ < (M + l)/2, and 



2 



■h 



N 
~2 



<^-Ia(m + q-i)-~(m-i), 

N 1 
<_-_A(M+Q-l). 



(13) 
(14) 



Here J s and Jh represent the uniform currents of down spins and holes, respectively. They 
are related with total momentum P tot = (2tt /N)(J S M + J^Q). The condition A=(odd 
integer) comes from the statistics of particles: down spins are hard-core bosons and holes 
are fermions. The three inequalities given above stem from the restrictions in mapping to the 
multicomponent Sutherland model. If the above conditions are satisfied, the wave function 
(|i~2|) yields the following eigenenergy: 

"a(a-i; 



N 2 1 1 
— E = -(A + 1) 2 M(M 2 - 1) - 2MJ S (N - J s ) + -(N 2 - 1) 

7T t 3 

+~(A + 1) 2 Q(Q 2 - 1) - 2QJ h (N - J h ) + 2\Q(J S - J h 



-N 
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~(A 2 -A-l)(iV 2 -l)Q 



+^q{a(M + Q)[(A + 1)M - AQ] + (A - 1) 



A.VQ-f (A + £) Q 2 -\ 



(15) 



When A = 1 i.e., Aj = 1 for all i, this reduces to the result by Wang et al. We obtain the 
ground state with magnetization by choosing both J s and as close as possible to N/2, 
while respecting the restrictions on M, Q, J s , and J^. 

Mapping (fj) to the multicomponent Sutherland model allows us to derive some excited 
states in addition to the states given here. Owing to the restrictions in the mapping, however, 
not all of states are found analytically in the anisotropic t-J model. The same applies to the 
long-range XXZ model. Thus the anisotropic models are different from completely soluble 
models: Bethe-ansatz solvable models or the 1/ sin 2 -type isotropic models. We emphasize 
that the anisotropic models belong to a new family of partially soluble models. 

We remark on the absolute ground state when the magnetization is varied. When A > 1, 
the state cannot be found analytically. Numerical calculation shows that the state has 
S* ot ^ with the parameter A larger than A c = 2.25 ± 0.05 for (N,Q) = (8,2). This is 
because the term with (A3 — 1) in (0) breaks the time reversal symmetry and acts like a 
magnetic field. In contrast to the spin model, presence of carriers in the t-J model should 
always lead to the gapless spectrum. We have also studied the case where the term with 
(A 3 — 1) in (g) is omitted. The absolute ground state has then S* ot = up to A = 10 for 
(N, Q) = (8, 2) and (8, 6) according to the numerical diagonalization. We remark that the 
ground-state eigenfunction is no longer of Jastrow-type without the (A3 — 1) term. 

We have seen that the Jastrow-type wave functions are exact solutions of the anisotropic 
models provided that the anisotropy parameters satisfy the solvable conditions. We now 
calculate the overlap | (^exactlV'j) I between the exact eigenfunction VWct and the Jastrow- 
type wave function ipj in order to see the validity of the latter function as a trial function 
beyond the condition for the exact solution. 

Figure |1] shows the overlap against the Ising anisotropy parameter A for 12 sites of the 
spin model. We first note that the overlap is almost unity over the whole range of the 
anisotropy. For example, its value is 0.97693 for A = 10 and Sf ot = 0. This means that 
the Jastrow-type wave function constitutes an excellent approximation even if the Ising 
anisotropy does not satisfy the condition for the analytic solution. In addition, the values of 
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A where the overlap is unity are: (a) A = 1; (b)A =1 and 6. These A satisfy the condition 
for the analytic solution. Second, comparing the top figure (a) with the bottom one (b), 
we find that the overlap in the case of Sf ot = 2 is closer to unity than that for 5 t z ot = 0. 
Finally, the overlap as a function of A is nonmonotonic with a conspicuous minimum and a 
faint second minimum. This behavior for S* ot = is common to other cases with different 
number of sites such as N = 6, 8 and 10. 

One may naturally ask whether the Jastrow-type function remains a good trial function 
if holes are introduced. Figure ^| shows the overlap against the anisotropy parameter A 
for the t-J model with 8 sites. Note that the anisotropy parameter is not Aj's but A. 
It is evident that the overlap is also very close to unity for the whole range of A. For 
A = 10, M = 1 and Q = 2, for instance, its value is 0.99965. In the case where A =1 and 
3 with (N,M,Q) = (8,1,2), the solvable conditions are satisfied. Hence the overlaps are 
unity in that case. In the limit of large A the t-J model describes localized holes and down 
spins. Then the ground state is a regular lattice of three components: up and down spins, 
and holes. As in the case of the spin model, the overlap shows a nonmonotonic dependence 
on the anisotropy parameter. 

One of the reasons why the Jastrow-type function is so good a trial function is the 
following: As A increases the wave function described by (|12|) becomes sharply peaked for 
evenly spaced coordinates of down spins and holes. The center of gravity is moving with 
a definite value of the total momentum. This state becomes the exact ground state of the 
Hamiltonian in the limit of large A. Here the ground-state energy does not depend on the 
total momentum, as seen by the absence of Jh and J s terms of the leading order 0(A 2 ) in 
(|15|). Thus the Jastrow-type function ( |12"D represents the exact eigenstate not only in the 



isotropic limit, but also in the limit of large A. It is still surprising that the Jastrow-type 
function interpolates the two limits so precisely. 

We now discuss the excitation spectrum of the present anisotropic model. One may 
ask whether the concept of spinons and holons remains useful in the case of anisotropic 
models. As the first step toward detailed investigation of the spectrum of the t-J model, 



we use numerical diagonalization to derive the excitation energy of the XXZ chain for each 
value of the total momentum K. Figure [| shows the spectrum for the system with 12 
sites and A = 10. It is evident that there is a quasi- continuum of excitations with energy 
gaps from the ground state. This continuum has a strong resemblance to the result of the 
nearest-neighbor spin model by Ishimura and Shiba ||. The almost vanishing excitation 
gap (= 3.3008 x 10~ 3 J) at K = tt reflects the degenerate ground state in the Ising limit. 

In summary, we have constructed the anisotropic t-J model which is analytically solvable 
for a restricted set of states. The model is probably not a completely integrable model unlike 
the other isotropic long-range models. We have shown that the Jastrow-type wave function 
is appropriate as the trial function over the whole range of the anisotropy parameter and 
filling. The excitation spectrum of the t-J model will be discussed elsewhere. 
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FIGURES 

FIG. 1. The overlap vs the Ising anisotropy A for the anisotropic spin model. The total 
azimuthal spin S^ ot = N/2 — M is: (a)5f ot = and (b)S* ot = 2. N and M denote the numbers of 
sites and down spins, respectively. The minimum and the maximum are shown by up and down 
triangle, respectively. 



FIG. 2. The overlap vs the anisotropy parameter A for the anisotropic t-J model. The total 
azimuthal spin S* ot = (N — Q)/2 — M is 2. N, M and Q denote the numbers of sites, down spins 
and holes, respectively. 



FIG. 3. The excitation spectrum of the anisotropic spin chain model with N = 12, A = 10 and 
S* ot = 0. E g represents the ground-state energy. 
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